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^— V ' We report on the transition between two regimes from several-atom clusters to much larger 

^vj , nanoparticles in Ar magnetron sputter deposition of WSi2, and the effect of nanoparticles on the 

properties of amorphous thin films and multilayers. Sputter deposition of thin films is monitored 
PfJ' by in situ x-ray scattering, including x-ray reflectivity and grazing incidence small angle x-ray scat- 

^ I tering. The results show an abrupt transition at an Ar background pressure Pc; the transition is 

associated with the threshold for energetic particle thermalization, which is known to scale as the 
product of the Ar pressure and the working distance between the magnetron source and the sub- 
strate surface. Below Pc smooth fllms are produced, while above Pc roughness increases abruptly, 
consistent with a model in which particles aggregate in the deposition flux before reaching the growth 
surface. The results from WSi2 films are correlated with in situ measurement of stress in WSi2/Si 
^^ ' multilayers, which exhibits a corresponding transition from compressive to tensile stress at Pc- The 

tensile stress is attributed to coalescence of nanoparticles and the elimination of nano-voids. 

PACS numbers: 61.05.cf, 68.35. Ct, 68.55.A-, 68.60.Bs , 68.65.Ac, 81.10.Bk, 81.15.Aa, 81.15.Cd, 82.30.Nr 

C3 ■ I. INTRODUCTION 

—^ ^ Magnetron sputtering has been widely used since the 1970s as a deposition method for metal, semiconductor, and 

^ . other inorganic thin films for applications such as optical coatings and microelectronic circuitsji Recently, there has 

Q ' been renewed interest in applications that require strict control over thin film structure and mechanical properties. 

O . For example, multilayers fabricated by sputter deposition are in development for x-ray optics such as in multilayer 

Laue lenses for nanometer scale focusing of x-rays3^— They require sub-nanometer roughness over several thousand 

/«f-\ , layers with minimal built-in stress, and these specifications have proven to be difficult to achieve with state of the 

^ ' art deposition techniques. Varying the background gas pressure can produce significant effects in thin films and 

OO multilayers, such as an abrupt change in interface roughness, or a sudden transition in film stress.i'^ For example, 

^D FuUerton et al7 have found that interfacial roughness in Nb/Si multilayers increases dramatically when the Ar pressure 

r ■ ' exceeds 9 mTorr. Cyrille et al.^ have observed a similar effect in Fe/Cr multilayers, and have made use of the pressure- 

^ ; induced roughness to enhance the magetoresistance in these structures. Similarly, striking changes in film stress have 

lO been observed for a number of materials. Hoffman and ThorntoniS have found that a stress transition from compressive 

^^ ' to tensile occurs at a pressure that depends on the atomic mass of the deposited material, and that the transition 

^— ^ ' pressure is inversely proportional to the atomic mass of the sputtering gasJ^ It would be beneficial to understand 

. . \ the origin of these effects because it would potentially lead to improved processes for fabrication of films for a variety 

^ • of applications, and because there is significant interest in understanding the fundamental mechanisms that govern 

surface dynamics during film deposition. 

In this paper, we demonstrate that the sudden changes in roughness and stress as the background gas pressure is 
C^ ■ varied both arise from an effect that is intrinsic to the magnetron sputtering process. There is a transition in the 
sputtered flux in which the dominant species produced by the magnetron source abruptly changes from atomic-size 
species to nanoparticles containing several hundred atoms. Our analysis methods involve in situ film deposition studies 
using synchrotron x-ray scattering to observe cluster distributions on the substrate surface before film coalescence 
and the power spectrum of surface roughness of films deposited from atoms vs from nanoparticles. The spectra are 
used to extract parameters related to the roughening processes on the film surface using standard equations of x-ray 
scatteringji^ The results show that surface roughening of deposited films at different pressures increases dramatically 
above a critical background pressure Pc- Analysis of the data shows that the increase in roughness cannot be explained 
by surface relaxation processes alone. Rather, it is caused by a sharp increase in the deposited particle volume, which 
results in increased roughness because larger particles contribute more deposition noiseJ^^Ai Specifically, we have 
investigated amorphous WSi2 film deposition using in situ synchrotron x-ray scattering and infer that there is a 
transition from atoms to nanoparticles in magnetron sputtering when the sputtering pressure is raised above Pc = 6 
mTorr. These results are described in detail in Sec. IIII Bl The data analysis are described in Sec. |IVl 
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We also report a complementary set of experiments involving measurement of thin film stress during the deposition 
of WSi2/Si multilayers. The results confirm that the stress transition occurs aX Pc — 6 mTorr, in good agreement 
with the results of the film roughness experiments. The tensile stress observed above the transition is interpreted 
as arising from coalescence of deposited clusters. The effect is similar to reports in the literature of tensile stress 
arising from coalescence of hills formed on the growth surface of amorphous films in the later stages of growth due 
to continuous viscous coalescence, and driven by surface tension due to the large curvature at the cusps between the 
hills. -'^^ This idea becomes more compelling with the realization that the deposition flux above the stress transition is 
fundamentally composed of nanoclusters, which may significantly contribute to tensile stress as they coalesce into a 
continuous film through the elimination of nano-voids. These results are presented in Sec. El 

This observation of a particle volume transition opens up possibilities for fundamental studies of nanoparticle 
aggregation, and applications such as film deposition and crystal growth via nanoparticle assembly. Production of 
clusters containing hundreds of atoms by aggregation of vapor has previously been demonstrated only in specialized 
gas aggregation instruments operating in a higher pressure range not typically used for sputter depositionJ^iii Our 
observations are for one specific material, however it is clear that the effect extends to a variety of other materialsi^ 
and possibly also to other deposition techniques such as pulsed laser deposition (PLD). In this regard, we note that 
PLD has some relevant similarities to magnetron sputtering, particularly the production of a dense vapor of energetic 
particles, which might produce clusters through the mechanisms of vapor-phase aggregation that we discuss in this 
paper. Furthermore, linking the stress transition and roughening transition to the particle volume suggests new 
possibilities for tailoring structural and mechanical properties of thin films and multilayers. 

II. EXPERIMENTAL DETAILS 
A. WSi2 thin fllm deposition and in situ x-ray scattering 

The growth experiments were performed in a custom-built ultrahigh vacuum chamber with a base pressure of 10"^*^ 
Torr installed at X2I station of National Synchrotron Light Source (NSLS), Brookhaven National Laboratory (BNL). 
A schematic of the experiment is shown in Fig. [IJa). We take the z direction to be sample normal and the j/-direction 
to be along the projection of the incident x-ray beams onto the sample surface. The x direction is vertical and 
the y-z scattering plane is horizontal in the laboratory frame. All amorphous WSi2 samples were prepared at room 
temperature by a dual-gun dc magnetron sputtering system (2" Meivac MAK) using ultrahigh purity (99.999%) Ar 
gas. A water-cooled 2-inch-diameter WSi2 target (purity 99.999%) was mounted on the gun, which was 90 mm away 
from the sample surface. The target normal lies in x-z plane and 7° to the substrate surface normal (z direction). 
Before deposition the target was pre-sputtered for 5 min with the target shutter closed. The dc sputtering power was 
kept constant at 50 W and the Ar pressure was adjusted ranging from 3 to 18 mTorr. The pressure was monitored 
by Pirani gauge (Stanford Research System), which was calibrated for Ar gas. The pressure varied by less than 0.1 
mTorr throughout the deposition process. Under these conditions, the resulting deposition rate v was between 0.103 
and 0.111 nm/s at all pressures. Amorphous thermal Si02/Si with rms roughness cr„ns = 0.25 nm were used as 
substrates for x-ray scattering measurements. 

The synchrotron x-ray flux at the X21 beamline is approximately 1 x 10^^ photons/s at photon energy of 10 keV 
(wavelength A = 0.124 nm~^) after the Si(lll) double-crystal monochromator and toroidal mirror. A slit of dimension 
0.2x0.5 mm^ (vertical x horizontal) 1070 mm in front of the sample and a slit (50x 1 mm^) 790 mm behind the sample 
define the angle of incidence a^ and the exit angle a/, respectively. The linear detector is oriented parallel to the 
sample surface. It consists of 384 pixels with 8 pixels/mm along the x direction. It covers a range of in-plane angle 
■0 = 2.88° at a distance of 955 mm behind the sample. The ultimate resolution of the instrument within the scattering 
plane is given by Sqz ~ 0.09 nm"-'^ for the direction perpendicular and Sqy ~ 10^'^ x q^ nni^^ parallel to the surface, 
while the resolution perpendicular to the scattering plane dqx ~ 0.007 nm^^. 

Four types of scans were performed during or after each deposition: (i) real-time grazing incidence small-angle 
x-ray scattering (GISAXS) monitoring of the evolution of surface morphology. The scattered intensity was measured 
by a one-dimensional detector oriented parallel to the sample surface. The exit angle a/ was kept constant at 0.6° 
and the incident angle was fixed at a^ = 0.2° to the sample surface. This asymmetric scattering geometry avoids the 
saturation of the linear detector because the specular beam does not reach the detector. In addition, the angle of 
incidence is below the WSi2 critical angle for total external reflection at 10 keV photon energy in order to enhance 
the surface sensitivity, and to keep g^ as low as possible. Each spectrum thus represents the scattered intensity as a 
function of q-x on a range of momentum transfer 0.003 nm~^ < Qk < 2.5 nm~^ at constant qz — 0.7 nm~^. This gives 
access to the surface roughness on lateral length scales between a few nanometers and 1 /im. (ii) In situ specular 
reflectivity scans after each deposited layer. This was done by rotating both the detector and sample about the x 
axis to keep the angle a^ = a/. The intensity was recorded by the sum of central 8 pixels of the linear detector. 
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FIG. 1: (a) Schematic of real-time grazing small-angle x-ray scattering (GISAXS) measurements during sputter deposition. 
The image is shown rotated with respect to the actual experiment, where the x direction is vertical. The incident x-ray beam 
(wave vector ki) impinges on the sample surface under a grazing incidence ai. The scattered intensity is recorded as a function 
of the out-of-plane angle a/ and the in-plane angle ij}. The wave vector transfer is denoted by Q = fc/ — fei. A 2D GISAXS 
reciprocal-space map is shown, which is acquired by varying a/ through the range 0.6°-10° while keeping ai = 0.2°. This scan 
is repeated at two detector positions to cover ^l> angle ranging from —0.2° to 5.56°. The negative q^ part is mirror imaged. 
This ctf scan closely approximates a q^ vs q^ map of the diffuse scattering component since qy is very small in this GISAXS 
geometry. For simplicity of notation we drop the subscript on q^ in reporting the results in the main text since the surfaces 
are assumed to be isotropic. The illustration also shows a TEM image of WSi2 film of nominal thickness 1 nm for a visual 
representation the sample surface in real space, [(b)-(e)] TEM images of a series of WSi2 films deposited on thin carbon grids 
at 8 mTorr at room temperature. The nominal film thicknesses are 0.5 nm, 1.0 nm, 2.0 nm, and 3.0 nm, respectively. The 
small white bar at the lower left of each image indicates a length scale of 2 nm. 

The surface roughness and thickness of these films were derived from those scans using a least-squares fit in which 
the reflected intensity was computed based on recursive application of the Fresnel equations;^^ (iii) in situ ip scans 
after each deposited layer. This scan circumvents the limitation of the length of the linear detector, which determines 
the range of Qx accessible in GISAXS with a fixed detector position. The linear detector is rotated about the sample 
normal by exactly the length of the detector at each step of the scan. For this study, we have employed scans to 
q^ « 10 nm~^. This extended qx range is particularly useful for quantitative modeling of the roughness spectrum, 
(iv) Two-dimensional GISAXS ip - af maps acquired while keeping a; fixed. This a/ scan closely approximates a g^ 
vs qx reciprocal map of the diffuse scattering component since qy is very small in the GISAXS geometry. This type 
of scan is shown in Fig. [TJa). 

B. Transmission electron microscopy (TEM) study of WSi2 cluster coalescence 

Amorphous ultrathin carbon coated copper grids (PELCO, No. 01824 from Ted Pella, Inc.) were used as substrates 
for TEM observations. TEM images were obtained on a JOEL 2100F high-resolution analytical transmission electron 
microscope operating at 200 kV at the Center of Functional Nanomaterials (CFN) at BNL. WSi2 clusters deposited on 
these carbon surfaces and observed via TEM were found to be in good agreement with cluster distributions observed 
via x-ray scattering on Si02 surfaces. Both carbon and Si02 are relatively inert, and therefore would not be expected 
to produce dramatically different results for our study, where surface diffusion and particle migration play a minor 
role. 

C. WSi2/Si Multilayer Experiments 

A well known technique to measure stresses in thin films is to measure the curvature and then to apply Stoney's 
equatioufi^ which in its biaxial form is given by^Sii^ 



1 /I- t^\ tf 



(1) 



Here the curvature is 1/i?, where R is the radius of curvature of a wafer, tf and tg are the thickness of the film and 
substrate, respectively, a is the stress in the film, E is the Young's modulus of the substrate, and v is the Poisson's 
ratio of the substrate. This equation applies to biaxial bending in the thin film limit, that is, in the limit in which 
the effective elastic constants are those of the substrate. We note that if the stress in the film is constant, then the 
prefactor in Stoney's equation is constant, and the above equation yields a linear relationship between curvature and 
film thickness. 

For thin films grown on Si substrate wafers, one must account for the fact that the elastic constants of Si are 
anisotropic in evaluating E and v.^^ However, for bending of Si(lOO) wafers, the biaxial modulus, E/{1 — ly), is 
conveniently isotropic in the plane of the wafer i^ 

Curvatures in WSi2/Si multilayers were measured in situ in a growth chamber located at Argonne National Labora- 
tory. Substrates were positioned on central vertical rotary axis, and sputtering proceeded horizontally. The chamber 
was equipped with a multibeam optical sensor (MOS) supplied by kSAi^^ The MOS was operated with a laser and 
beam splitters to produce a 4x3 (columnxrow) grid of beams reflected from a multilayer sample. The sample curva- 
ture was measured in both the horizontal and vertical directions. The MOS was situated on a vacuum port adjacent to 
the sputtering gun, and curvature measurements were made by periodically bringing the sample from a position facing 
the sputtering gun to a position facing the MOS. This was done at intervals corresponding to 1.1 nm of thickness. 

The substrates were 50 mm diameter Si(lOO) wafers. Multilayers consisting of 20 bilayers of WSi2/Si with layer 
thicknesses of 5.5 nm each. Multilayers were studied for the Ar plasma pressures of 2.3, 6, 12 and 18 mTorr. These 
parameters were chosen to be the same as for an earlier in situ x-ray reflectivity study^Sl and layer thicknesses were 
checked subsequently with ex-situ x-ray reflectivity measurements made after a full multilayer was grown. 

III. RESULTS OF WSI2 THIN FILM DEPOSITION EXPERIMENTS 
A. Transient stage of WSi2 deposition 

Figure[Ha) displays a 2D GISAXS map of WSi2 clusters deposited on thermal Si02/Si substrate at 8 mTorr for 10 
s, which has a nominal fllm thickness of 1 nm. Symmetrical rounded shapes of the scattered intensity are observed, 
indicating the presence of three-dimensional clusters on the growing surface. The peak of the intensity occurs at 
Qpeak ~ 2 nm~^, indicating a particle separation of 2Tr /qp^ak ~ 3 nm, and the width of the peaks indicates that the 
cluster positions are nearly randomly distributed. 

Figure [Ijb)-(e) shows TEM images of WSi2 clusters on amorphous carbon grids with different nominal fllm thick- 
nesses from 0.5 to 3 nm. They were deposited at 8 mTorr with deposition rate of 0.105 nm/s. Clusters 1-2 nm in 
diameter are seen at 5 s, which corresponds to about 0.5 nm average thickness. At later times the clusters aggregate 
into elongated meandering islands on the surface, indicating limited coalescence. Similar clusters are also observed at 
4 mTorr (not shown in the figure). The TEM observations are consistent with the GISAXS pattern. 

These observations clearly show that clusters are present on the surface during the early stage of film deposition 
before coalescence. One of the key questions that we endeavor to answer in this study is whether the clusters 
form in the sputtering plasma before reaching the substrate, or whether the aggregation occurs on the surface. Both 
possibilities are plausible. Aggregation in the gas phase depends on collisions of energetic particles with gas molecules, 
which produces a thermalization effect. This process is expected to be pressure dependent (for reference, we note 
that the mean- free path in Ar at 8 mTorr is about 6.8 mm). On the other hand, surface aggregation depends on 
the mobility of particles after they land on the surface. Processes such as thermal diffusion should be very limited 
when the substrate is held at room temperature, and so the length scale of aggregation in this case is expected to be 
short. More generally, both gas aggregation and surface aggregation processes may act in concert to produce the final 
structures observed. We believe that this is the most likely explanation of the fact that we observe clusters on the 
surface at all pressures above and below Pc- Therefore, the observed cluster size of 1 - 2 nm observed at the earliest 
stage of deposition should be taken as an upper limit of the cluster size deposited from the magnetron source. We 
will return to the question of deposited particle size and surface cluster size in the next section, and also in the data 
analysis presented in Sec. IIVI 

B. Pressure dependence and temporal evolution of surface roughness 

The rms roughness derived from specular x-ray refiectivity is plotted in Fig. [5] as a function of Ar pressure. The 
films deposited at different Ar pressures each have the same thickness tf — 180 nm. The surface roughness shows a 
very abrupt transition at 6 mTorr where the roughness increases dramatically. It is well described by a power law 
(P - PcY, with Pc = 6 mTorr and s = 0.30 ± 0.05. 



E 
c 



o 
cr 
w 



EC 1 



Ot 



— 1 1 1 r- 



— 1 1 1 r- 




_i I I i_ 



_i I I i_ 



5 10 15 

Ar Pressure (mTorr) 



20 



FIG. 2: The rms roughness of f/ = ISO-nm-thick WSi2 films as a function of Ar background pressure in which two roughness 
regimes are observed. The red fine is a power law (P — Pc)" ■ 



A striking dependence of surface morphology of thin films and multilayers on the sputtering gas pressure has been 
observed previouslyj^i^ i^^'^^ At low gas pressures, sputtered and back-reflected species ballistically impinge on the 
growing surface with higher average kinetic energies, which leads to enhanced surface relaxation and smooth surfaces. 
Surface smoothing by energetic particle bombardment at low pressures is consistent with our data, as we discuss below; 
see Fig. |4]and the discussion of mechanisms in Sec. IVI Al The number of gas-phase collisions in the path from target 
to substrate increases with pressure so that above a "thermalization pressure", the velocities of energetic particle 
emerging from the source are reduced to thermal velocities before they strike the substrate j ^'^'i^^ Our suggestion for 
explaining the roughening transition, which we discuss further in Sec. IVI A( is to combine this thermalization effect 
with a model of cluster growth in the gas phase. This model predicts a very sharp onset of cluster growth and hence 
roughening as the sputtering pressure is varied, in agreement with the results shown in Fig. [51 

Before proceeding, we mention several alternative models for the roughening transition based on changes in the 
strength of smoothening processes with constant particle size: (i) a possible model is that the surface roughens above 
the thermalization pressure simply because energetic smoothing processes are suppressed above Pc- This model seems 
highly plausible. However, this model is inconsistent with the shape of the curve in Fig. [5] In particular, a gradual 
variation in the strength of the smoothing processes as a function of the background Ar pressure is predicted, which 
would produce a continuous increase in the roughness. In contrast, we observe a sudden onset at Pc- (ii) A second 
possibility is a transition from stable surface kinetics at low pressures to an unstable surface at higher pressures, 
which can be understood in terms of a change in sign of one of the smoothing coefficients in a linear continuum 
model.— A sharp transition near Pc is predicted as a surface instability mechanism becomes dominant.—^— In this 
model, the surface would transition to an unstable surface that would be characterized by rapidly increasing amplitude 
at a preferred spatial wavelength. In contrast, as we describe below, the surface is found to exhibit stable kinetic 
roughening in both pressure regimes below and above Pc- Therefore, alternative models for the roughening transition 
based on changes in the strength of smoothening processes alone with no variation in particle size are inconsistent 
with the data presented. 

Figure [3] shows real-time GISAXS evolution during the deposition at 8 and 4 mTorr, respectively. In both cases, 
with growing amount of WSi2, a transient peak associated with 3D clusters first appears and moves towards lower q, 
indicating coalescence. This is consistent with the TEM results in Fig. [IJb) . Film coalescence is complete at about 4 
nm in each case, which begins the kinetic roughening stage of the surface evolution. Below, we discuss Fig. EJa) and 
|3Jb) in turn. 

Figure [3ja) illustrates the typical evolution of the GISAXS intensity for pressures above Pc- After the film coales- 
cence the diffuse scattering profile S{q) increases monotonically with film thickness for q < Qc, where qc is the cutoff 
wave number. As detailed in Sec. IIV Al after a correction at low q, the GISAXS data are proportional to the power 
spectral density (PSD) which is expected to be constant for q < qc-, i.e., the spectrum is expected to be white below 
a cutoff wave number. We note that qc is related to the surface correlation length ^, the maximum length scale to 
which the surface roughness has propagated, by qc — 27r/^. The intensity increases (at least) linearly with time since 
S{q) cc fl tf "''^ for q < qc, where fi is the particle volume and tf is the film thickness, and typical values for the 
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FIG. 3: Real-time GISAXS evolution during sputter deposition (a) at 8 mTorr and(b) at 4 mTorr. The arrows indicate the 
transient peaks observed in the precoalescence stage. 



exponent (2 + 2a)/z range from 1.0 for linear models such as the Edwards- Wilkinson (EW) model^^ to ^ 1.75 for 
the Kardar-Parisi-Zhang model<^ See Sec. IIVI for further discussion of the scattering intensities. We interpret this 
rapid increase in the intensity to be indicative of a large particle volume. The spectrum also exhibits characteristics 
of kinetic roughening because the curves for different film thickness overlap perfectly at g > gc )^^'^^ which also rules 
out the possibility that the surface is unstable. At later times, a crossover between scaling regimes exhibiting different 
exponents is observed at ~ 1 nm~^. We interpret this as being due to two distinct relaxation processes for length 
scales larger and smaller than this characteristic length because of the change in the slope of the curve above and 
below 1 nm~^. Note that this feature does not shift after the coalescence, which suggests that it is the result of 
steady-state processes rather than being simply a remnant of the transient feature. 

Figure injb) illustrates the behavior below Pc- After coalescence, the film surface is statistically similar to that of 
the starting surface. This is because the shot noise is very low due to the smaller particle volume, making the increase 
for q < Qc and the cutoff itself very difficult to observe at low sputtering pressures. At later times the spectrum is 
nearly identical to the one shown for tf — A nm even for much thicker films. The q~^ dependence of the curve is the 
stable spectrum, indicating an EW-type smoothing mechanism. The roughening transition is therefore interpreted as 
a transition from an EW regime at low pressures to a new kinetic roughening regime above Pc that is triggered by 
the larger particle volume. 
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FIG. 4: (a) GISAXS spectra show the PSD profiles after deposition of the smoothing layer at 5 mTorr (red) on rough layer 
deposited at 8 mTorr (blue). The substrate spectrum is also shown for comparison. During the deposition of the low-pressure 
layer, the roughness decreases from 1.6 nm to 0.66 nm. The smaller wavelength surface roughness decays away faster and the 
corresponding roughness spectrum exhibits a power law with an exponent of —2; (b) GISAXS intensity evolution at various q 
in a semilog plot. The high spatial frequency features vanish faster. The slope of the curves shows the decay rate, which is 
plotted in (c). (c)Plot of the decay rates, confirming that the relaxation is dominated by a q'^ term, as indicative of the EW 
model. The a2 relaxation coefficient is also extracted from this curve. 



In Fig. m an experiment to determine the smoothing coefBcients is shown, which consists of depositing a relatively 
rough 100 nm film at 8 mTorr, and then smoothing the surface by depositing 220 nm at 5 mTorr. This data can 
be used to extract surface relaxation coefficients directly from the smoothing kinetics. We have previously used this 
technique to study relaxation mechanisms in ion bombardment induced smoothing of surface.-22. The corresponding 
GISAXS profiles for Si02 substrate, rough and smoothed WSi2 layers are shown in Fig. SJ^a). It can be seen that 
the diffuse intensity decreases during deposition of WSi2 at 5 mTorr, consistent with the decrease in the rms surface 
roughness, which was determined by specular x-ray reflectivity. Note that the increasing intensity for q < O.l nm"-'^ is 
due to the fact that the initial GISAXS curve for the rough layer does not accurately give the spectrum of roughness 
because the small roughness condition is not satisfied due to the large initial surface roughness (see Sec. lIVCp . The 
spectrum of the smooth layer exhibits a dropoff at qc, the cutoff frequency, which is related to the characteristic 
correlation length on the surface corresponding to the largest surface features which have already been smoothed out. 
Above this spatial frequency, a scaling behavior with an exponent of -2 is observed. This part of the spectrum is 
identical to the one obtained for WSi2 deposited on smooth Si02 at pressure below the transition at Pc, as shown in 
Fig. ^h). This behavior confirms that the stable roughness spectrum is independent of the starting surface. 

Figure nib) shows the same GISAXS intensity plotted vs time at various length scales. The curves follow a decaying 
exponential form, in good agreement with the prediction of linear theory (Sec. IIV A[) . The roughness decay rate at 
different wavelength is extracted by an exponential decay fit for all the curves. The decay rates are plotted as a 
function of spatial frequency in Fig. \M,c). It shows unambiguously that a g^-dependent smoothing mechanism is 
dominant on the surface during the deposition below the transition pressure. Thus, EW-type behavior is confirmed 
for pressure below Pc- 



10 F 



10 



10 



10' 



10' 



10 



10 



I I I I I 1 1 ll I I I I I 1 1 ll I II 




60 nm WSij @ 4 mTorr 

1 nm WSij© 18 mTorr 

2 nm 
4 nm 

1 nm ° »»'™— 

30 nm „„„_ . 

50 nm 
lOOnm 
180nm 

1 I I I ■ ■ ■ 



0.01 0.1 1 

q ( nm"^ ) 



FIG. 5: Real-time GISAXS evolution during sputter deposition at 18 mTorr on a smooth layer deposited at 4 mTorr. 
arrow indicates the transient peak associated with 3D clusters on the film surface. 
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Figure [5] shows a different variation in the experiment. Here the first layer is deposited at 4 mTorr, producing a 
very smooth WSi2 surface, then a second layer is deposited at 18 mTorr. There is no driving force for clustering on 
the surface when depositing WSi2 onto itself, such as a difference in surface energiesi^i^^ Therefore, by observing 
the real-time GISAXS spectrum we can study the transition from smooth to rough without the complication factor 
of depositing onto a surface that promotes clustering. The results show that the transient peak associated with 
nanoclusters on the surface appears again at ~ 2 nm~^. Note that the transient peak is not part of the kinetic 
roughening, but rather is due to the deposition of isolated clusters on the substrate. This is the stage of deposition 
before the particles begin to interact, and was discussed by Edwards and Wilkinsoni^ This observation clearly confirms 
that the clusters observed at P > Pc are generated by the sputtering source rather than by aggregation on the surface. 

The bump at 1 nm~^ for tf > 30 nm in the steady-state spectrum may be related to a surface instability process. 
However, the surface still exhibits kinetic roughening because smoothing processes dominate over the instability during 
the deposition. The discussion of roughness dynamics is in Sec. IIVI 



IV. THEORY BACKGROUND AND ANALYSIS OF ROUGHENING DYNAMICS 

A. Surface morphology dynamics 

Combining the surface relaxation mechanisms with the stochastic roughening within a hnear approximation, we 
can write a kinetic rate equation of the Langevin type, for surface growth: 

^h^^-h{q)-h{q,t) + ri{q,t), (2) 

4 

6(g)=«^a„.g", (3) 

where /i((7,i) is the Fourier transform of surface height h(r,t), q = {qx,qy) is the in-plane wave number and g = ^n = \q\ 
is the circular averaged wave number. The factor v is the deposition rate. The term h{q) on the right of Eq. ^ 
models the surface relaxation processes during film growth. Coefficient a„ is a constant characteristic of the specific 
lateral mass transport mechanism indicated by n. The term ri{q,t) is the reciprocal-space stochastic noise term that 
describes the random arrival of the depositing species. It has the property 

{7iq,t))=0 , {r,{q,tUq',t'))=vnSiq + q')6{t-t') , (4) 

which represents uncorrelated deposition noise with strength vfl. The factor v is the deposition rate, and Q is the 
volume of the species being deposited, whether they are atoms, molecules, clusters, or nanoparticles. 

The PSD of a growing surface is defined as (|/i(g,i)p), which is also the Fourier transform of the correlation function 
C{r, t) = {h{r, t)h{r' , t)). The radially averaged PSD can be analytically determined from Eq. ^, 

PSD{q,t)^{\h{q,t){')^PSD{q,0)e-''('^^* + vn^^^^-^ , (5) 

PSD{q,t)^[PSD{q,0)^PSD{q,oo)]e-^'''^'i^' + PSD{q,oo) , (6) 

where PSD{q, 0) is the power spectral density of the substrate, which equals to zero for perfectly smooth surface. 
The first term on the right of Eq. ([S]) damps the contribution from the substrate, while the second term presents 
the roughness which increases from intrinsic growth processes. Stochastic roughening by random deposition creates a 
surface that contains features of all sizes. The PSD is just a constant fine for q < q^, i.e., the q spectrum is white. The 
constant value depends on the rms roughness and increases as fl tf for any linear model of the type represented by 
Eq. ([2]). ri is the particle volume in deposition fiux and tf is the film thickness. Eq. ([7]) shows in the long-time limit 
(t — > oo), the spectrum reaches the steady state that decreases as a power law g~", or multiple power laws if there is 
more than one nonzero a„, depending explicitly on the identity of the relaxation mechanisms. In the framework of 
kinetic roughening, the surface morphology exhibits asymptotic scaling behavior at long times and long length scales. 
In the case of Family- Vicsek scaling, 92, the PSD can be expressed by the relationi^ 

PSD{q,t)^q-^''^+^^f{q/q,) , q^^T^/^ (8) 

with the scaling function 



/(") 



{const, if u 3> 1, 
u2"+2 ifu<l. 



where a is called the roughness exponent, z the dynamic exponent, and /? — a/z the growth exponent. The integral 
of PSD{q) over q allows the time dependence of the surface roughness, 

/■oo 
aL.W- / dq^PSD{q,t) , (9) 



which yields a power law arms ~ t^ ■ Therefore, the PSD evolution contains the complete spatial and temporal 
descriptions of the surface morphology evolution. Note that for EW surface dynamics, a = 0, /3 — 0, z = 2, and so: 



At) 



87ra2 



-In 



(10) 



where tr, 



P /a2 is suggested by Natterman and Tang, and I is a minimum length scale 
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B. Structure factor of a rough surface 



Sinha et al. have showr^ that within the distorted- wave Born approximation the structure factor for an isotropic 
surface is given by 



Siq'z,q\\) 



^Re(q'^fal 



dr (el^^l'^M - l) Jo(g||r) , 



(11) 



where q'^ = k\ \/v?- — cos o?^ + \ln^ — cos o?, J denotes the vertical momentum transfer in the sample, fc — 2tt / \ is 



the wave number, n is the index of refraction, and cFrms is the rms roughness, as determined from specular reflectivity 
data as shown in Fig. [21 C{r) is the surface correlation function with C(0) = <^1ms^ ^^^ ^o denotes the Bessel function 
of the first kind and zero order. 

In the case of small roughness when (|g^|crrms)^ ^ 1, the integral of Eq. (|lll) can be approximated by a Fourier 
transform of C{r), which is the PSD of a rough surface. Therefore, the real-time x-ray scattering gives a direct way to 
monitor the surface morphology evolution during the film growth over the entire spectrum of accessible length scales. 



C. Extraction of the particle volume and relaxation coefficients from S{q) 



(a) Rough surface S(q) and PSD 
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(b) Smoothed surface S(q) and PSD 
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{c) PSDs extracted from S(q) data 




FIG. 6: Extracted PSD from normalized GISAXS intensity profile for film deposition at (a) 8 mTorr on thermal Si02/Si 
substrate; (b) subsequent deposition at 5 mTorr on as-grown rough surface; (c)the extracted PSD evolution for the smoothing 
experiment, which is consistent with linear theory, i.e., see Eq. ((5)1. 



When {\q'^\arms)'^ ^ 1? ^^^ small roughness approximation is not valid because we can not make the approximation 
gkjl C{r) _ ]^ ~ |q^pC(r) in Eq. ([TT|). Then, the x-ray scattering intensity does not accurately gives the PSD of 
a rough surface. In particular, we find that the scattered intensity underestimates the PSD at the lower values of 
q. However, the PSD can be recovered by following procedure: (i) numerically reverse Fourier transform the x-ray 
intensity in Eq. pT|) to get the function F(r) = e'^-' '-^'^''^ — 1; (ii) normalize F(r) according to F(r)|r=o = e'''^' '^r,ns — l^ 
since C(0) — o^rms'^ (iii) calculate the correlation function from C{r) = ln{F{r) -{- l)/\q'J^'^] (iv) numerically Fourier 
transform C{r) to get the corrected power spectral density. 

In Fig. [HJa), circles represent the x-ray scattering intensity normalized with roughness arms for a rough surface 
deposited at 8 mTorr, which is also shown in Fig. HJa). The red line in Fig. IH^a) is a smooth curve obtained according 
to Eq. ([SJ to fit the experimental data. It is used to numerically calculate the corrected PSD, which is shown in 
squares. The green line is the fitting of PSD to extract the smoothing coefficients, as well as the particle volume in 



10 



TABLE I: Film parameters, particle volumes and smoothing coefficients. 



Sample 


p 


tf 


Cfrms 


n 


02 


a-i 


a4 


\/a2/a4 




(mTorr) 


(nm) 


(nm) 


(nm«) 


(nm) 


(nm^) 


(nm^) 


(nm-i) 


I (A) 


8.0 


100 


1.60 


3.55 


0.34 


-0.25 


0.42 


0.90 


I (B)" 


5.0 


220 


0.66 


0.038 


0.10 


0.00 


0.00 


— 


11(A) 


8.0 


480 


2.84 


3.55 


0.12 


-0.25 


0.42 


0.53 


II (B)" 


4.0 


600 


0.62 


0.020 


0.14 


0.00 


0.00 


— 


III 


5.9 


200 


0.39 


0.038 


0.12 


-0.08 


0.08 


1.19 


IV 


6.5 


188 


1.52 


1.80 


0.25 


-0.05 


0.85 


0.54 


V 


7.0 


183 


1.80 


2.50 


0.28 


-0.10 


0.80 


0.59 


VI 


8.0 


181 


2.10 


3.60 


0.35 


-0.40 


0.85 


0.64 


VII 


10.0 


182 


2.57 


6.62 


0.50 


-0.50 


1.00 


0.50 


VIII 


14.0 


179 


3.15 


N.A.' 


N.A. 


N.A. 


N.A. 


— 



"Subsequent deposition atop layer (A). 
''Integration loses accuracy when qzCFrm 



> 2. 



the deposition flux. Fig. ^h) shows the same data analysis for the subsequent smooth layer growth at 5 mTorr. Note 
that the calculation of the PSD for the overgrowth at 5 mTorr (layer B) atop the 8 mTorr layer (layer A) requires 
both sets of relaxation coefficients and particle volumes for layers A and B, i.e., the final curve includes the history 
of the first layer. The two corrected PSD curves for the rough surface at 8 mTorr and smooth surface at 5 mTorr are 
also plotted in Fig. |6l[c) for direct comparison. The data correction and fitting procedure were repeated for several 
films deposited at different pressures, and the extracted smoothing coefficients and particle volume are listed in Table 
III 




6 8 10 

Pressure ( mTorr ) 



FIG. 7: Particle volume produced by sputtering source extracted from GISAXS spectra. The line is a power law (P 
for the transition with an exponent of m. = 0.70 ± 0.10 and transition pressure at 6 mTorr. 



Pc) 



The particle volume, H. from Table U is displayed vs sputtering pressure in Fig. [71 Clearly, there are two main 
regimes separated by a sharp transition at ~ 6 mTorr. The particle volume changes by two orders of magnitude as the 
transition boundary is crossed. A single WSi2 molecule has a volume of 0.042 nm"^, so the particle volume above the 
transition pressure 6 mTorr corresponds to several hundred atoms, while below 6 mTorr the deposition flux mainly 
consists of single atoms and/or very small (few atoms) clusters. The particle volume is described by a power law 
(P - PcY" with exponent of m = 0.7 ± 0.10. In Fig. H we have shown that arm.s -- (P - PcY with s = 0.30 ± 0.05, 
so that m « 2s. Comparing to Eq. (jlOp . we see that (7rms oc vs2 is expected for linear theories of kinetic roughening. 
Evidently, arms is dominated by the change in particle volume, and is only affected slightly by changes in the strength 
of relaxation processes. 
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FIG. 8: Curvature data plotted as a function of sputtering time for Ar plasma pressures of 2.3 mTorr (solid squares), 6 mTorr 
(open circles), 12 mTorr (solid circles), and 18 mTorr (open squares). Twenty bilayer periods of thickness 11.0 nm each were 
deposited at each pressure. The raw data were shifted vertically to agree at the growth start. 
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FIG. 9: Curvature values extracted from Fig. [8] at the conclusion of growth, that is, at 10000 s at which time 220 nm had been 
deposited, as a function of sputtering pressure. 



V. STRESS EVOLUTION RESULTS FOR WSI2/SI MULTILAYERS 



in situ curvature data plotted as a function of sputtering time for different Ar pressure are shown in Fig. [S] 
As indicated in the figure, compression produces a negative curvature (convex deformation) according to our sign 
convention for curvature, and vice- versa for tension. The compressive film stress develops at low sputtering pressures, 
while tensile stress is observed for high sputtering pressures. The curvature values at the conclusion of growth are 
plotted as a function of pressure in Fig. |9l It clearly shows a transition from compressive to tensile stress with 
increasing sputtering pressure. The transition pressure is at 6 mTorr, in good agreement with the results of the film 
roughness experiment. It suggests that the transition from compressive stress to tensile stress is correlated with the 
transition from smooth to rough film growth. 

The stress transition is very similar to those reviewed in Fig. 1.39 in the text by Freund and Suresh for the metals Cr, 
Mo, Ni, and Ta, and they are consistent with models in which tensile stress has been linked to island coalescence»iii^ 
The compressive stress observed below the transition can be interpreted as a result of " atomic peening" , which is due 
to the energetic particles striking the growing film with their high impact kinetic energy. 
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VI. DISCUSSION OF MECHANISMS 



This work combines several topics that have developed independently: the kinetics of surface evolution and dynamic 
scaling theory, cluster condensation in the gas phase, and pressure-dependent roughness and stress transitions of films 
during magnetron sputter deposition. In this section, we will discuss the mechanisms for each process. 



A. Plasma-gas condensation 
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FIG. 10: Model for pressure dependence of cluster growth incorporating the effects of thermalization and aggregation (line). 
The parameters used are Pc = 4.8 niTorr, Vd/no = 8.5 x 10~^° rn^/s, and Nc.o = 1. With these parameters, the cluster size in 
the high pressure limit, P — s- oo in Eq. (|13|l . is found to be 833. The data points are the cluster volumes from Table [J divided 
by the volume of a single WSi2 cluster. 

Current models for inert gas aggregation of nanometer scale clusters postulate a nucleation step followed by cluster 
growth due to atomic vapor condensation onto the newly-formed embryos i^^i^^ In order for nucleation to occur, 
the kinetic energy of the sputtered atoms must be reduced to thermal energies via collisions with the background 
gas. While plasma-gas condensation instruments are typically operated at high pressures of 200 - 300 mTorr where 
thermalization occurs over distances of a few millimeters^ in low-pressure sputter deposition thermalization occurs 
over a distance that can be comparable to the source - sample distance. Somekh^^ has performed Monte Carlo studies 
of this process and has found that the thermalization can be characterized by a pressure-distance (PD) product that 
gives the critical pressure at a given drift length. Their results predict that W sputtered atoms with an initial kinetic 
energy of 5 - 25 eV lose 90% of their energy with a PD in the range of ^ 120 Pa-mm (900 mTorr mm). At P = Pc, 
PD= PcL, where L is the distance between the sputtering source and the sample. Given our drift length of L = 90 
mm, a critical pressure of 10 mTorr is predicted, which is in reasonably close agreement with our observation for WSi2 
of Pc = 6 mTorr. Therefore, our results are consistent with the point of view that the thermalization of sputtered 
vapor is a controlling factor in determining Pc for cluster growth, as well as the roughness and stress transitions. 

The picture above implies that particles leaving the sputtering source travel ballistically towards the substrate 
surface until they are thermalized, and only then begin to coalesce into clusters. At this stage, the density of the 
sputtered vapor must equal a density no so that the clusters grow to a size, Nc atoms per cluster. Hihara and 
Sumiyama^ have developed a model based on collision cross sections in the vapor phase to describe the growth of 
clusters in plasma-gas condensation instruments at high pressures. We will adapt this model for low-pressure sputter 
deposition in order to approximately relate the sputtered metal vapor density uq to the particle size. Integrating the 
expression given by Hihara and Sumiyama for dNc/dz, where the z axis is in the direction from the sputtering source 
to the sample, neglecting re-evaporation from clusters, we have 



no = 



■nrl VthL 



(12) 



where vth = {^ksT /'KmaY''^ is the thermal velocity of the metal atoms in the vapor, r^ is their atomic radius, and 
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nia is the atomic mass. Also, Vd is the drift velocity of the clusters. Note that in our calculations we take a WSi2 
subcluster to be the fundamental unit rather than individual atoms for simplicity. 

From Table U we see that the particle size for 10 mTorr deposition is about 6.6 nm^^. This corresponds to Nc = 
156 WSi2 subclusters. The thermal velocity of a WSi2 subcluster is 163 m/s, and r^ = 0.216 nm, uia — 3.99 x 10~^^ 
kg. Using the thermal velocity of a fully developed cluster with Nc — 156 , we have as the drift velocity Vd = 13.0 
m/s, which yields an estimate of uq « 9.7 x 10^^ m^"^. This is a remarkably high vapor density given that the density 
of the background Ar at 10 niTorr and 300 K is 3.2 x 10^° m"'^. As a check of the consistency of this result with our 
experiment, we also relate this number with an estimate for the flux of sputtered particles. Given / = 100 mA and 
V = 500 V on the target, a target diameter of 2 in., and assuming a sputtering yield of approximately one WSi2 per 
incident ion, we estimate J ^ 3.1 x 10^" m~^-s~^. Combining this with the relation J = n^Vd and substituting into Eq. 
(H^, we find that rig ~ 4.8 x 10"'^^ m^"^, which is close to the a priori estimate obtained from Eq. p^ alone. Therefore, 
although there is significant uncertainty in the parameter Vd on the order of a factor of 2, these estimates suggest 
that a sputtered metal vapor density of 10^^-10^° m~^ is required for production of clusters in sputter deposition. 
Since the corresponding telltale effects of tensile stress and roughness have been reported frequently in the published 
literature on sputter deposition of metal films, a very high sputtered vapor density must be inherent to the magnetron 
sputtering process. 

Equation (J12p does not predict any significant pressure dependence. However, it is possible to construct a model 
combining the effects of thermalization and cluster growth by substituting an effective drift length L^ff for the total 
drift length L where ie// is the net drift length over which clusters form, i.e., ie// = L — Lth- Here, Ltu is the 
thermalization length Lth — PcL/ P^ so that L^ff — L{1 — Pc/P) ■ We note that ie// should be zero at P = Pc- 
Rearranging Eq. (I12[) . substituting Leff for L, and explicitly including the embryo cluster size Nc_o, we have 



^c«<^ r/ , - .,. ^., P.^ , ,.i/3l3 . „. „ (13) 



Nc.o 




forP < Pc , 


{irrlnQVthL/Svd) (l - 


P ) ^ ^^c,Q 


for P> Pc 



This model successfully predicts a strong pressure dependence near Pc, as we show in Fig. [TO] The curve matches 
the experimental values well, except for within ~ 1 mTorr of the transition. This relatively good agreement suggests 
that the model captures the essential mechanisms responsible for the cluster transition. We emphasize that although 
a model can be constructed based on the pressure dependence of particle impact smoothing without any cluster 
aggregation, such a model will not resemble our particle volume vs pressure data because the smoothing processes 
vary continuously with particle energy, i.e., they do not exhibit a threshold at a definite pressure. Finally, we comment 
that the exponent used in the fit of Fig. [7] is not likely to be a universal exponent. Rather, the detailed shape of the 
cluster size vs pressure curve appears to be the result of at least two processes operating in concert, and thus may 
change under different experimental conditions. 

B. Kinetic roughening 

As discussed in Sec. IIV A) kinetic roughening is the result of the competition between the inherent noise in the 
growth process due to the nonuniform nature of the incoming fiux and surface relaxation effects. Each smoothing 
mechanism affects the morphology differently and leaves a different signature in the PSD of the growing surface. 

The roughening in the low-pressure regime is exemplified by the GISAXS spectra shown in Fig. |3l^b), and the 
smoothing layer deposited on a rough first layer shown in Fig. Hl^a). Strong smoothing is observed, and both curves 
exhibit a crossover to stable q"^ behavior above a cutoff wave number qc that shifts toward lower q with time. The 
fact that this behavior is independent of the starting surface is consistent with kinetic roughening, as can be seen in 
Eqs. (O - ([T]). The classic linear model involving only the 02 coefficient in Eqs. ([2]) and dS)), corresponding to a 
V^/i term in real space, was originally described by Edwards and Wilkinson for the case of particle sedimentationi^ 
EW behavior has been observed before in sputter deposition at low pressures4^^ More recently, this term has been 
suggested as being related to energetic impact induced downhill relaxation, which is expected to occur when energetic 
particles reach the growth surface with energies exceeding ^20 eV j^^i^^ 

The observed GISAXS spectra are more complicated in the higher pressure range. However, the detailed fitting 
results give some insight into the mechanisms. The 02 term listed in Table U is related to the EW (V^/i) term, which 
results in the q~^ behavior at q > qc- At first sight, it is somewhat curious that the 02 coefficient increases significantly 
for pressures above Pc because we would expect that impact-induced smoothing would be significantly suppressed 
above the thermalization pressure. However, relaxation of clusters into nearby hollows in the surface is also described 
by an EW-type term. Thus, we interpret the increase in the 02 coefficient as the appearance of a new relaxation 
process related to cluster relaxation on the growth surface. The monotonic increase in a2 with pressure is attributed 
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to the increasing particle size. The 04 term can be interpreted similarly as a shorter length scale process related to 
the sintering of clusters. For example, processes such as viscous flow can give a V^/i term4i The 04 term can also be 
related to MuUins surface diffusion,— although surface diffusion may be negligible at room temperature. The 03 term 
has a negative value, which shows that it is not a relaxation process. Rather it could indicate a surface instability 
process such as strain-driven diffusion. ^^ However, we emphasize that we have not found evidence for unstable surfaces 
under any of the experimental conditions, indicating that b{q) > for all q. 

C. Stress transition 

The stress transition is a particularly striking effect that has been demonstrated for a variety of sputter-deposited 
films.— 1"^°'^^ In poly crystalline films tensile stress observed at pressures above the transition has been suggested to be 
due to closing the small gaps between islands to form grain boundaries, which ultimately merge into a continuous 
fllm^^i^ 

A similar concept has been suggested as a way to explain tensile stress in amorphous films, which can form cellular 
structures characterized by grooves at the boundaries between hillocks J^s^'-^ The growth of the cellular structures 
is driven by the deepening of grooves due to shadowing effects, which can contribute to tensile strain by closing of 
the grooves. However, we point out that the coarsening behavior observed for these layers is distinct from kinetic 
roughening because larger features increase in size at the expense of smaller features and hence there is no steady- 
state power spectrum over any range of length scalesi^ Therefore, this type of coarsening is only relevant to a later 
stage of the film deposition, while our observations are for the early stage when shadowing effects are absent and the 
groove stress mechanism does not play a dominant role in producing stress. The fact that the stress transition in 
our data is coincident with the formation of clusters strongly suggests that the clusters directly produce the stress 
by sintering and elimination of voids between particles, which can be very fast for nanometer scale particles even at 
room temperature.^^ This shrinkage effect is well known in powder metallurgy and is due to that fact that as particles 
fuse together their centers also move closer together. ^■^ 

To conclude this discussion, we point out that the observation of clusters in magnetron sputtering suggests many 
opportunities for tailoring the properties of thin films. Aside from the control of the film growth process by varying the 
background pressure, other parameters, particularly the sputtering power and the type of sputtering gas are expected 
to have a strong effect on cluster formation and in the properties of films thus produced. 

VII. CONCLUSIONS 

In conclusion, we have developed a comprehensive model that encompasses several effects occurring in sputter- 
deposited thin films and multilayers. Specifically, we have shown that the roughening transition and the stress 
transition both arise from aggregation in the sputtering plasma as pressure increases above a threshold value. 
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